Three-body calculation of the rate of reaction p + p + e — > d + v e in the Sun 
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Using expansion of the three-body wave function of the pep system in the initial state on hy- 
perharmonic functions the rate of the p + p + e~ — > d + v a reaction in the Sun is calculated. The 
results of calculation of the flux at 1 AU are compared with the results of measurement made by 
the Borexino collaboration and Bahcall et al. theoretical predictions. 
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I. INTRODUCTION 

The study of the neutrino from the Sun is a multi- 
purpose problem. Indeed, the measurement of the solar 
neutrino flux incident on the Earth helps to clarify the 
properties of the neutrino, for example, the phenomenon 
of the oscillations and to determine its parameters (the 
mixing angle and the eigenmass of neutrino). Then, the 
value of the flux and its spectral properties contain in 
some cases information on the nuclear reactions that in 
the visible future cannot be observed in the laboratory. 
We mean here the reactions: 



and 



P + P 



P + P 



d + e+ + v e 



(1) 



(2) 



Moreover, as emphasized by John N. Bahcall in his book 
and in a series the papers [3-11] (see also @) solar neutri- 
nos bring the information about processes in the center 
of star connected with solar structure models. 

Although the reaction ^j, called the pep reaction, 
plays no significant role in hydrogen burning of stars, 
it is esential in detecting monoenergetic neutrinos of 
E v = 1.442 MeV. Also, a measurement of the neutrino 
flux from reaction ([2]) can be useful for the determination 
of the parameters of the Standard Solar Model. 

All above considerations contain sets of parameters 
used for fitting observable data in the framework of dif- 
ferent models like star models, nuclear reactions models 
and so on. An example of such a type is the two-body 
model [7[ for reaction (2) which is essentially three-body 

1- 

The purpose of this work is to exclude the model ele- 
ments as much as possible in the description of reaction 
([2]) and treat the initial state as a purely three-body state. 
From a strict point of view, the modern treatment of any 
nuclear reaction inside the Debye-sphere is still the model 
in a sense of the absence of the dynamical consideration 



(it is the six-body problem in the Sun interior condition) 
of particles inside the sphere. We hope, however, that the 
three-body instead of the two-body treatment of reaction 
([2]) is a step in the right direction. 

Below, we concentrate ourselves mainly on the process 
([2]) for the following reasons. Firstly, in 2012 the results 
of the first experimental observation of process ((2|) was 
announced @ , after more than 50 years of studies of solar 
neutrino problems. The second reason is connected with 
the absence of the three-body treatment of the initial 
state in the process ©■ 

There is a question concerning sensitivity to the choice 
of the nucleon-nucleon potential and related to these po- 
tentials wave functions of the bound state of the deuteron 
and the continuum state of the pep system. 

With all this in mind we present below the treatment 
which takes into account all peculiarities of the pep three- 
body system. In Section [TTJ we start from the inputs 
for considered problem: the weak Hamiltonian and NN- 
potentials, then in Section HlIl we consider solution of the 
Schrodinger equation to deremine the pep wave function 
of the initial state. In Section IIV1 we present the calcu- 
lation of the probability of the pep reaction and the as- 
trophysical S pep factor taking into account the Coulomb 
and strong interactions simultaneously. The results for 
the rate of the process and fluxes of the neutrino arc 
discussed in Section fVl and conclusions are presented in 
Section ED 



II. INPUTS 

The electron capture by the nuclear system can be de- 
scribed by the following non-relativistic effective weak 
Hamiltonian nil 



Hi, 



1 



7T 



-(+) 



l-g'Ut ^ (-) 



V2 



G v l • U 



i=i 



'Electronic address: irgaziev@yahoo.com 



G A cr ■ tji - Gpcr ■ uxCTi ■ u x ]8(r — n), (3) 

where i/\ = vjv {v is the neutrino momentum); 1,1$, er 
and <Tj are the 2x2 matrix unit operators and spin angu- 
lar momentum operators for the lepton and ith nucleon; 
r and ri are the space coordinates of the lepton and zth 



nucleoli; t^ + \ t\ ' are the isobaric-spin operators which 
transfom a lepton electron state into a lepton neutrino 
state and ith nucleon proton state into the ith nucleon 
neutron state; Gy, G a and Gp are the vector, axial vec- 
tor and "induced" pseudoscalar coupling constants, re- 
spectively. We take G v /{hcf = 1.153 x 10- n GeV" 2 
and G A /(hc) 3 = -1.454 x Kr n GeV~ 2 [11]. We can 
simplify the weak Hamiltonian for the pep system: the 
last term in Eq. (j3|) can be neglected because the emit- 
ted neutrino has energy E v = 1.442 MeV and this term 
encloses factor v/2m p , where m p is the proton mass; the 
pep — > d + u e transition satisfies the Gamow- Teller selec- 
tion rule; therefore, the first term of Eq. Q does not give 
contribution to the matrix element of transition; and at 
last we take into account that the electron neutrino has 
spin opposite to its momentum v. Finally, for the weak 
Hamiltonian we get 

A 

H w = t^G a 2 rH<7 • cn5(r - n ). (4) 

i=l 

With this Hamiltonian of Eq. Q we obtain the electron 
capture transition matrix element. 

The energy of thermalization of particles in the interior 
of the Sun corresponds to E ~ 1.3 keV which is small 
on the nuclear energy scale. Therefore, we can use the 
simple AW-potential which describes correctly the low 
energy data of the nucleon-nucleon system. We apply 
the Gauss and the Yukawa potentials [l2|]. We fit the 
parameters of these potentials to get the values of the 
deuteron energy, the scattering lengths and the effective 
ranges for pp and pn scattering. 

For the Gauss potential 

V(r) = -V Q exp(-r 2 /R 2 N ). (5) 

the calculation with the fitted parameters 

V S = 30.36 MeV, R S N = 1.816 fm 

gives the scattering length s a pp = —7.884 fm and the 
effective range s r pp — 2.678 fm for pp scattering at the 
singlet state (s), while the parameters 

V * = 60.572 MeV, R* N = 1.65 fm, 

lead to the scattering length t a pp — 5.484 fm, the effective 
range t pp = 1.85 fm for np scattering at the triplet state 
(i) and the binding deuteron energy Ed — 2.225 MeV. 
The second potential is the Yukawa potential 

V(r) = —^-exp(-r/R. N ) (6) 

with the parameters for the singlet state: 

V S = 44.05 MeV, R S N = 1.206 fm, 
and for the triplet state: 

= 53.27 MeV, R% = 2.43 fm. 



The results of calculation of the scattering lengths and 
effective ranges are: 

s a pp = -7.782 fm, s r pp = 2.868 fm, 

l a np = 5.626 fm, V np = 1.895 fm. (7) 

To find the flux of neutrino we must use some Standard 
Solar Model (SSM). There are several SSM which are in 
good agreement with the helioseismologically determined 
sound speed, temperature and density of elements as a 
function of solar radius, the depth of the convective zone, 
the surface helium abundance, and so on. We applied 
data of parameters presented in the model BS05(OP) 
[llij . The results of Bahcall et al. [13] shows that the 
flux from the pep reaction is not sensitive to the type of 
SSM. 



III. THE WAVE FUNCTION OF THE pep 
INITIAL STATE 

We note that Bahcall and May Q used the pep wave 
function in a factorized form as the product of the wave 
function of the relative motion of two protons and the 
wave function of the electron moving in the Coulomb 
field of these protons. However, such a representation 
is not a quite correct procedure due to the long-range 
nature of the Coulomb interaction, even for the asymp- 
totic behavior of the wave function when electron is at 
a large distance from the protons [HI, [l6| ■ All the more, 
we cannot perform such a factorization at a small rela- 
tive distance where we need to know the wave function 
with a sufficient accuracy to get a good accuracy of the 
calculation for the transition matrix element of the pro- 
cess pep — > dv e . Also, there is a problem with the to- 
tal angular momentum of pep, because the moments of 
subsystems are not conserved. Fortunately, the relative 
angular moment of two nucleons in the initial and final 
states is zero as angular moments of the electron and the 
neutrino; therefore, the last problem did not arise in the 
Bahcall calculations of the pep reaction. 

We use the hyperspherical harmonics [lj], [3 for ex- 
pansion of the pep wave function in the initial state and 
solve directly the three-body Schrodinger equation and, 
therefore, we are free of the problems which we mentioned 
above. As in the nonrelativistic approach the orbital and 
spin moments are conserved independently of one an- 
other, we can expand the spatial part of the pep three- 
particle wave function over hyperspherical functions and 
we obtain the linked system of the radial differential 
equations. Derivation of the system radial equations is 
given in Appendix |XJ To proceed further, we now use 
the system of Eqs (|A14[) for the radial wave functions. 
Since the total energy of the pep system is low, the main 
contribution to the three-particle wave function gives the 
states with the zero relative orbital moments. Owing to 
the centrifugal potential and up <C 1 at small distance 
where we need to calculate the wave function with higher 
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accuracy, contributions of the components with the hy- 
permoments K > should be suppressed in the total 
wave function. Taking into account these conditions, we 
need to find the solution of the single equation only for 
K = l x = l y = with omitting the non-diagonal terms. 
We omit all the indices because they correspond to zero 
values of the quantum numbers and get the equation for 
the radial wave function U(p): 



d 2 U(p) 



+ 1 -'f-(v(p) + ±- K 2 )u(p)=0, (8) 



1 dU 

p dp v • p- 
2p23,E > (E is the total energy of the pep 



dp 2 

where k 2 
system); 

V(p) = V N (p)+V c (p), 
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Here 773 is the three-body Coulomb parameter which is 
defined as 
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The matrix element of the nuclear potential is the 
following: for the Gauss potential 
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and for the Yukawa potential 

16M23Vb 3T^- / 2(i 7 )+i 2 (^ 
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N 



where I n (z) is the modified Bessel function of the first 
kind, and L n (z) is the modified Struve function. 

To find a unique solution of Eq. ([8]) for the continu- 
ous state of the pep system, we must determine bound- 
ary conditions. Instead, to define the function U{p) and 
its derivative at the origin (p — 0), we define the wave 
function in a point p close to zero because we know be- 
haviour of the wave function near the origin. We can 
take boundary conditions as 

U{p ) = J2(k po), U'(po) = kqJ^oPo), (17) 



where k = \/ n 2 + \ V N (po) I- 

The wave function U(p) at large distance, where the 
nuclear interaction is negligible, has the following asymp- 
totics: 

U(p) e iS * cos 5 3 (F 00 (k P ) - tan <$ 3 GooM), (18) 
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FIG. 1: (Color online) The pep radial wave function (solid 
curve) obtained by solving Eq. (J5J and its asymptotics at 
p S> Rn (dashed curve). The pure Coulomb wave is dot- 



dashed curve. E pep = 6 keV. 



where 63 is the three-body scattering phase shift, 
Fqo(kp) (Goo(kp)) is a three-body regular (irregular) 
Coulomb wave function. The Coulomb wave functions 
have asymptotics [13 at the large distance 

/ 2 / 5 \ 

Fqo(z) = \ — cosfz - 773 In 2z - -it + a c ), (19) 
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where z = up and a c — argT(5/2 + %r\ 3 ) is the three- 
body Coulomb scattering phase shift when hypermoment 
K and the total angular moment L are equal to zero, and 
T(z) is the gamma function. 

We solve Eq. ([8]) using the boundary conditions (TT7)) 
and then matching the logarithmic derivative of the solu- 
tion in the asymptotic region with the logarithmic deriva- 
tive of the asymptotic solution (Eq. (|T8"|0 we define the 
three-body phase shift £3 which depends on the total en- 
ergy E of the pep system. Of course, this phase shift is 
very small because npN <§; 1, where pn is equal to the 
radius of NN interaction. Matching the condition of the 
numerically obtained solution of the Schrddinger equa- 
tion with its asymptotics also allows us to find the nor- 
malization factor. At great distances the ratio of the un- 
normalized solution to the asymptotic function becomes 
a constant. 

Figure [T] shows the results of calculation of the pep ra- 
dial wave function (solid curve) at E = 6 keV and its 
asymptotics for the Gauss NN potential (dashed curve). 
Since the value of the total energy of the system is small, 
the scattering phase will be small too. Therefore the 
asymptotics of the solution can be described by a pure 
three-body Coulomb wave (dot-dashed curve). We ob- 
tain the same results for the Yukawa NN potential. We 
need to know the three-body wave function with high ac- 
curacy in the interval < p < 35 fm because the deuteron 
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wave function cuts integration at the distance 35 fm when 
we calculate the matrix element of the pep — > d + v tran- 
sition. 



IV. THE PROBABILITY OF THE pep^d + v 
REACTION AND ASTROPHYSICAL S FACTOR 

Using the weak Hamiltonian (Eq. Q) we can write 
the matrix element of the pep reaction as 

(+) 



H if =G A <Vv \ (TT (+> I if e > X 
2 

< I I %ep >, 



(21) 



where ip e and ip v are the spin functions of the electron 
and neutrino, respectively; vf^ is the wave function of the 
deuteron, while ^ pep is the three-body wave function of 
the pep system in a continuous state. 

Since the weak interaction is small, the first-order per- 
turbation theory can be applied for calculation of the 
probability of reaction per unit of time P3; therefore, we 
have 



2tt- 



P* = -r-\ H„ 



~ 2 P{E„) 



(22) 



where the overline means summation over spins in the 
final state and averaging over spins in the initial state, 
and the density of neutrino states is 



P(E U ) 



2tt 2 H 3 c 3 ' 



(23) 



Here E v , c are the neutrino energy and light speed, re- 
spectively. Using the energy conservation law we define 



E = E + E 



(M d - 2m p - m e )c 2 



where E pp is the kinetic energy of proton-proton relative 
motion, Md, m p , m e are the deuteron, proton and elec- 
tron masses, respectively. If we put kinetic energies of 
the particles to zero, we get E v = 1.442 MeV. 

Making averaging in Eq. (|22|) and using the Jacobi 
coordinates for the particles defined by Eq. (|A5[) of Ap- 
pendix |A] we obtain 



Pa = 



3E?,G 



/ tf3(xi)¥pep(a:i,yi)| 



(24) 

Note, the presence of the delta-function in the weak 
Hamiltonian ((4]) leads to calculation of the integral over 
variable y\ to the significance of the pep wave function 

at the distance t/io — 



mi(m 2 +m 3 ) /r, ^ / ro e 

M23 (mi +m 2 +m 3 ) x ' ~ \J 2m/ 1 ' 

The deuteron wave function can be written in the form 



- ^Voo(*i), 



•X'l 



(25) 



where the radial wave function must have a correct value 
of the asymptotic normalization coefficient to get a cor- 
rect value of the probability of the reaction. Using Eq. 
(|25p and taking into account only the first term in the 
expansion of the pep wave function over the hyperhar- 
monics functions we obtain for the overlap integral 



16^ 



(26) 



Here the three-body radial wave function is calculated at 
the point p = \J x\ + y\ ~ x\, and ^3 is the three-body 
nuclear scattering phase shift modified by the Coulomb 
interactions. 

Usually, in the nuclear astrophysics for parameteriza- 
tions of a two-body reaction cross section the astrophys- 
ical S factor is used [l9j]. To find the S factor from the 
cross section at a small energy region, the Gamow factor 
is extracted from the cross section, i.e. 



(27) 



where r\ = Zi^e 2 /(hv) is the Sommerfeld (Coulomb) 
factor for the colliding nuclei with charges Z\ and Zi 
having the relative velocity v. A similar procedure can 
be carried out for the reaction with three particles in the 
initial state because the three-body radial wave function 
Ukl (p) m tne continuum state at a distance larger than 
nuclear interaction radius encloses factor exp(— 77773/2) | 
T(K + 5/2 + 1773) |. Therefore we define the astrophysical 
S pep factor for the pep reaction as 

P 3 (E) = G (E)S pep (E), (28) 



G (E) 
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(29) 



S pep (E) is almost a linearly varying function of E and 
S pep {0) is not equal to zero, like the astrophysical S pp 
factor for two-body reactions. Note, the unit of S pep 
coincides with the unit of J3 because Gq(E) is dimen- 
sionlcss and it is the Gamow factor for the pep reaction. 
If we use expansion of S pep (E) over E as 

S pep (E) = S + S 1 E + S 2 E 2 , (30) 

we obtain the following results for the value of the coef- 
ficients: 

For the Gauss potentials 

50 - 2.38 x 10 10 fm 6 /s; 

51 = 3.03 x 10 10 fm 6 /(MeVs); 

^2 = 1.45 x 10 10 fm 6 /(MeV 2 s); (31) 

and for the Yukawa potentials 

50 = 2.33 x 10 10 fm 6 /s; 

51 = 3.01 x 10 10 fm 6 /(MeVs); 

5 2 = 1.78 x 10 10 fm 6 /(MeV 2 s). (32) 
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(a rate constant) is 



ICm-.n — < P3 >— 



-pep 
00 00 oc 



<p(v e ) tp(v pi )(p(v P2 )P 3 (E)dv e dv Pl dv P2 , (36) 



000 



where <p(v) = ^{^f 2 exp(-f^) is the Maxwell- 
Boltzmann distribution function for particles of mass m 
and velocity v and k is the Boltzmann constant. Exclud- 
ing the velocity of the center of mass of the pep system 
by transferring to the relative velocities and the fact that 
the center of mass of the pep system is almost the same 
as the center of pp we get 



FIG. 2: (Color online) The pep astrophysical S pep factor as 
a function of energy calculated by the Gauss (solid line) and 
Yukawa (dashed line) NN potentials. 



K--pep — 



Lp(v e )(p(vpp)P 3 (E)dv e dv P p, (37) 



Behavior of the S pep astrophysical factor on energy is pre- 
sented in Fig. [2] where linear dependence of the S factor 
is seen clearly for both potentials. The small difference 
between the S pep factors for the Gauss and Yukawa po- 
tentials can be explained by small differences of the ob- 
tained low energy scattering parameters. To check the 
validity of the results of our calculations of the pep re- 
action, we calculated the astrophysical S pp factor for the 
pp reaction, according to the conventional definition (Eq. 
(|27[)) by the same NN potentials. If S pp is approximated 
by a polynomial 



is the 



where tp(v pp ) = 47rt$ p (^ r ) 3 / 2 exp(-^J , v pp 

relative velocity of protons and fj, pp — m p /2. Taking into 
account that 



y 2 
■ pp 



E = 



and making transformations 



Vpp — Vcosa, v e = .l^^-Vsma, E 



H PP V 2 



we obtain 



Spp(E) — So + S\E + S2E 2 



we get the following results: 
For the Gauss potentials 



5 = 4.21 * xlCT 25 MeVb; 

51 = 4.75 x l(T 24 b; 

5 2 = 3.02 x 10- 23 MeV _1 b; 



and for the Yukawa potentials 



50 = 4.11 x l(T 25 MeVb; 

5 1 = 4.64 x l(T 24 b; 

5 2 = 2.93 x 10" 23 MeV _1 b. 



(33) 



(34) 



(35) 



These results are very close to data presented in 19]. 



V. RATE OF THE pep REACTION AND THE 
SOLAR NEUTRINO FLUX 



The number of the pep reactions per unit of time per 
unit of the solar interior volumes of protons and electrons 



Kpep ~ (fcT) 3 



-E/kT 



P 3 (E)E 2 dE. (38) 



Using Eqs. (|28j). ([29]) we finally obtain 



v 

We note that the rate constant depends on the temper- 
ature, the nature of the reactants, but does not depend 
on their concentration. If we define the Gamow energy 
as 



E G - 



1 f2Tre 2 z e ff\ 2 



(40) 



2 V h 

the integrand of Eq. (|39[) has the maximum at the energy 



E, 



%ax = (^-kT^/Ea^ 



2/3 



(41) 



Here z eff is defined through Eqs. (fi"2"l). (JT3J) and JH]) and 
it equals to 



16 , 

z eff = ^-( a i + a 2 + 03). 
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(42) 
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FIG. 3: (Color online) The pep (solid line) and pp (dashed 
line) rates versus the radius of the interior of the Sun. The 
results are presented for the Gauss AW potentials. The be- 
haviour of the rate curves for the Yukawa AW potentials is 
very close to the displayed curves. The results are presented 



in units of 10 7 (pp) and 10 5 (pep) cm 
of the Sun. 



Rq is the radius 



Note, if we take z e f f = 1 we obtain the point of maximum 
of the integrand in the equation of the pp reaction rate 
constant. 

To obtain the rate of reactions, the rate constant must 
be multiplied by the density of the reactants, 



n 



pep 



(43) 



where n e and n p are the numbers of electrons and protons 
in the unit of volume. 

All solar parameters like temperature, densities of pro- 
tons and electrons vary with the radius of interior of the 
Sun. The results for the neutrino fluxes presented by 
Bahcall et al. [f| show that the fluxes from the pp and 
pep reaction are not sensitive to the considered solar mod- 
els. Therefore, for calculation of the rates of the pp and 
pep reaction we apply the BS05(OP) model @. 

Figure [3] shows the rate of the pp and pep reaction as 
a function or the solar interior radius for the Gauss AW 
potentials. We obtain the same behaviour in the case of 
the Yukawa AW potentials. We see that the rate of the 
pep reactions is more than a hundred times less than the 
rate of the pp reaction for the whole distance from the 
center of the Sun, and reactions occur in the core of the 
Sun where the temperature and density are higher. 

Integrating the reaction rate (|43|) over the volume of 
the Sun we find the total flux of neutrinos emitted by the 
Sun. Dividing this total flux by the area of the sphere of 
the radius of one astronomical unit (AU) we obtain the 
neutrino flux $ passing through a unit area of the Earth 
surface. The results of calculations of the neutrino fluxes 
are presented in TableQ] Note, the results for the neutrino 
flux from the pp reaction are close to the results obtained 
by Bahcall et al. while there are differences between the 
fluxes from the pep reactions. 



TABLE I: Predicted fluxes <E> PP and <& pep (without survival 
probability), in units of 10 10 (pp), 10 8 (pep) cm _2 s _1 . 



Standard Solar 
Model 

BS05(OP) 



BS05(OP) 



6.20 2.04 



6.05 1.99 



304 



304 
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our results 
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potential 
our results 
with Yukawa 
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1.45 


417 
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The Borexino collaboration announced the results of 
the neutrino flux measurement: <fr pep = (1.6 ±0.3) x 
10 8 cm _2 s _1 [9|. Taking into account that the survival 
probability of the neutrino (due to the neutrino oscilla- 
tion) in the pep reaction predicted by the Borexino col- 
laboration equals P = 0.62 ± 0.17 at 1.44 MeV [§|, we 
find that the neutrino flux at 1 AU should be equal to 



pc p = (1.27 ± 0.35) X 10 8 cm- 



for the Gauss 



potential, and * pep = (1.24 ± 0.34) X 10 8 cm~ 2 s _1 

for the Yukawa potential. 

We note that our results for the neutrino fluxes from 
the pep reaction obtained by taking into account the sur- 
vival probability lie within the confidence interval of the 
experimental data. At the same time, the Bahcall re- 
sults are out of this limit at all fluxes listed in Table 
UJ Comparing our calculated low-energy parameters for 
the Gauss and Yukawa potentials, we see that they differ 
by 2-7%, and the neutrino fluxes from the pp have 2% 
difference too, and the results of $ pp of Bahchall et al. 
shown in Table Q] for all SSM are not more than 2-4% 
from our calculated Therefore, we may conclude 

that the dependence of the results on the type of AW 
potentials is very weak. However, the difference in the 
results obtained by us and Bahcall et al. can reach up to 
39-45%, which means strong sensitivity to the choice of 
the wave functions of the initial three-body state of the 
pep system. 



VI. CONCLUSION 

In the framework of the three-body approach the prob- 
ability of the process pep — > d + v e under the condition 
of the solar core has been found. The rate of the above 
process and neutrino flux are found and compared with 
the Borexino experiment and previous calculations. The 
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value of the neutrino flux obtained from the pep reaction 
in the three-body treatment appeared to be ~ 40% larger 
as compared to the Bahcall et al. value. This can be un- 
derstood as a correct description of the movement of the 
electron producing the screening effect between protons. 
In this work we have introduced the astrophysical S fac- 
tor for the three-body reaction which is an analog of the 
S factor introduced for binary processes. To discriminate 
between different star models on the basis of our results, 
it is necessary to essentially reduce experimental errors 
in the above experiment. 
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coordinates the Hamiltonian Hq of the free motion of 
particles can be written as 

Ho = (A Xl +A m )= (A^ + A^ 

M23 v ' 1^23 K 



M23 



A* 3 + A y3 



(A6) 



We define the square of the hyperradius and hyperangle 
as: 

2 2, 2 2, 2 2, 2 ^rn \ 

Xi = pcosai, yi = psmcti, on G [0,tt/2]. 

In the variables of hyperradius and a set of angles f2^ = 
Hi, on) [xi and t/j are unit vectors determining the 
azimuthal and polar angles)), we can rewrite the operator 
Hq as 



Appendix A: Schrodinger equation 

The Hamiltonian of the pep system is 
H = H + V 23 + V 123 , 
where Ho is kinetic energy operator: 



Ho = —. 



2// 



-A 



23 



23 



Al(23), 



(Al) 



(A2) 



subscript 1 means electron, while 2 and 3 correspond to 
two protons; fi 23 is the reduced mass of two protons, 
while M 1 (23) is the reduced mass of the system of two 
protons and electron; 



V c — V C 4- V C 4- V c 
v 123 ~ V 12 ' v 23 ' v 31 



(A3) 



is the sum of the Coulomb potentials in the pep system, 
and is nuclear potential of interaction between two 
nucleons. 

The wave function of the initial state of the pep sys- 
tem is the eigenfunction of the Hamiltonian (IA1|) . which 
has a continuous spectrum of energy and satisfies the 
Schrodinger equation 



HV(ti,T2,t s ) = £#(ri,r 2 ,ra). 
Let us define the Jacobi coordinates: 



(A4) 



rrijmk 



-(Tj - r k ) 



(rrij + m k )n23 

m^rrij + m k ) 
(mi +m 2 + m 3 )p 23 
mjTj + m k r k 



—r. 



rrij + m k 



(A5) 



where fi 23 = m 2 m 3 /(m 2 + m 3 ) is the reduced mass of 
two protons, indices ijk=123, 231, or 312 and to, (rj) 
is the mass (coordinate) of particle i. In introducing 



tt _ n 2 ( d 2 5 d 1 2( 



where the operator iC 2 (f^J is 

K 2 (n,) = -_-4cot2a 4 — + 



(A7) 



1 



1 



cos z a 



, l 2 (x i ) + — 2 —l 2 (y i ), (A8) 



sm on 



where I 2 = — Ag iV is an angular part of the Laplace op- 
erator. 

The hypersphcrical function is defined as the solution 
of the equation [ljj 

K 2 <S>(a t ) = K{K + 4)<f>(Qi), K = 0,1,2,3,... (A9) 

The quantum number K is called hypermoment, the 
eigenfunction $(f2j) is 

(A10) 

where the function $^.» i « ma: * rn! '< (fij) is 

^{% m ^ myi m=N l ^ lv *(cosa i ) l ^(sma i ) l n x 

P> +1A ^ +1/2 (cos2ai)^ 4 m iCj (*i)^ jro , j (yi); (All) 

(ixiTnxJyiTnyi \ LMj is the Clebsch- Gordon coefficient; 
Yi^m^x) is the spherical function; 



N 



l Xi ly 

K 



I 2nl(K + 2)(n + l Xi + l yt + 1)! 
r(n + l x , + 3/2)r(n + l y . + 3/2) 



(A12) 



n = (l/2)(K—l Xi — l Vi ) must be an integer number; P n ! 
is Jacobi polynomial. 
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Let p and q be conjugate momenta to the coordinates 
x and y. Then determining square of the wave number 
K 2 = 2/123-E/^ 2 (E is the total energy of the pep system 
in the cm. frame) for the pep system in the continuous 
state, removing a free motion of the center of mass of the 
pep system from Eq. (|A4j) and using the following expan- 
sion of the pep wave function of a continuous spectrum 
over the hyperharmonics functions 



* P , q (x,y) = (27r) 



■K u KL 



(p) 



KLM \ 



M 2 



(0 P )$^(Q K ) (A13) 



we get the radial Schrbdinger equation: 



d 2 U l KL V (p) 
dp 2 



dp 



{K + 2f 



U 1 K 1 L V {P) 



E 



;1 kk'lm(p) u k'Up)> 



(A14) 



where superscript (*) denotes taking complex conju- 



gate of the function, and & K l M {Q P ) = ®klm(^i)> 
$^£ M (il K ) is defined on a hypersphere of unit radius in 
a six dimensional momentum space. The matrix element 



V 



KK'LM 



j(p) equals to 



V 



KK'LM 



}{P) 



2/x 



2.; 



h 2 



*l x ly 

KLM 



(v£(\xi I) + vg(\ Xl \) + vg(\x2\y 



(A15) 



i'X. 



\))®K>lM( n p) dn , 



The system of Eq. (|Al4j) is a system of the linked one 
dimensional equations which must satisfy the boundary 
conditions depending on the particular physical situa- 
tion. 

It is easy to show that near p 



regular solutions of 
(|A14|) must have the behaviour as p K+2 . If we omit the 
nondiagonal terms in the equations, solutions at p — > oo 
have the asymptotic behavior as a superposition of the 
regular and irregular Coulomb functions [17| . 



[1] J. N. Bahcall, Neutrino Astrophysics, Cambridge Univer- 
sity Press, 1989. 

[2] J. N. Bahcall, |arXiv: astro-ph /0209080] 

[3] J. N. Bahcall, and M. H. Pinsonneault, Phys. Rev. Lett. 
92, 121301 (2004). 

[4] J.N. Bahca ll, A.M. Serenelli, and S. Basu, ApJ, 621, 
L85 (2005); |arXiv:astro-ph/0412440] 

[5] J.N. Bahcall, A. M. Serenelli, and S. Basu, ApJ Suppl., 
165, 400 (2006). 

[6] |http://www.sns.ias.edu/~jnb/| 

[7] J. N. Bahcall and R. M. May, ApJ 155, 501 (1969). 

[8] V.B. Belyaev, S.V. Levin, S.L. Yakovlev, J. Phys. B37, 
1369 ( 2004). 

[9] G. Bellini et al. (Borexino Collaboration), Phys. Rev. 

Lett. 108, 051302 (2012). 
[10] H. Primakoff, Rev. Mod. Phys. 31, 802 (1959). 
[11] B. Povh, K. Rith, C. Scholz, F. Zetsche, Particles and 

Nuclei: An Introduction to the Physical Concepts, Sixth 



Edition, Springer, 2008. 
[12] J. M. Blatt, V. F. Weisskopf, Theoretical nuclear physics, 

New York: Springer- Verlag, 1979. 
[13] |http://www.sns.ias.edu/~jnb/SNdata/Export/BS2005/| 

bs05op.dat. 

[14] J.N. Bahcall, A.M. Serenelli, Sarbani Basu, ApJ 621, 

L85 (2005); ApJ Suppl. 165, 400 (2006). 
[15] E. O. Alt and A. M. Mukhamedzhanov, Phys. Rev. A 47, 

2004 (1993). 

[16] A.M. Mukhamedzhanov and M. Lieber, Phys. Rev. A 
54, 3078 (1996). 

[17] R. N. Djibuti and K. V. Shitikova, Method of hyperspher- 
ical functions in atomic and nuclear physics, Moscow: 
Energoatomizdat, 1993 (In Russian). 

[18] M. Fabre de la Ripelle, S. A. Sonanos, R. M. Adam, An- 
nals of Physics 316, 107 (2005). 

[19] C. Angulo, et al, Nucl. Phys. A 656, 3 (1999). 



8 



